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ie LNCRODUCT ION 


This investigation considers a non-linear technique for 
extracting natural resonances from transient electromagnetic 
scattering responses of radar targets. These natural reso- 
nances represent the complex poles of the target's transfer 
function in the Laplace transform s-plane. The advantage of 
their use in identification is their dependence only upon the 
geometry and composition of the target and not upon the aspect 
and polarization of the incident signal. 

Concepts and methodologies evolving from this idea have 
been developed based on the Singularity Expansion Method 
(SEM) by Baum [Ref. 1] in which the location of s-plane 
Singularities characterizes a target. The direct extraction 
of poles (and residues) from the complex impulse response was 
initially attempted by Main and Moffat [Ref. 2]. These two 
early references spawned the idea of using natural resonances 
as the basis of non-cooperative target recognition (NCTR). 

The application of the Householder orthogonalization 
method using a successive orthogonalization process and the 
eigenvalue method, that finds eigenvectors and eigenvalues, 
can determine the number and value of poles from the tradi- 
tional transient signal model of a scatterer if there is a 
sufficiently low level of noise in the data [Ref. 3]. But 


both methods show poor quality of performance if we apply 








them to the new signal model, even were there is no noise 
assumed in the signal data. This new signal model, which 
represents the electromagnetic scattering transient response 
of actual targets, was recently derived by Morgan [Ref. 4]. 
As considered in [Ref. 2], the relationship between a 
scattering target and the singularities being extracted from 
the natural complex resonance signal has been emphasized in 
the context of the number of poles and the accuracy of poles 
to ensure the unique relationship between an object and its 
image. Namely, the accuracy of pole extraction plays an 
important part in setting baselines to check consistency 
between a received signal and a set of singularities. Also, 
[Ref. 2] describes the modeling problem of a signal such that 
we express a transient response signal in terms of complex 
exponentials and the way it can be manipulated by Prony's 
algorithm. This signal model (termed here as the traditional 
model), was found not to represent the actual transient signal 
Bimeenemrecent work by Davenport [Ref. 5]. In 1982, Manilha, 
in his work, commented that the new signal model might be 
Pamehetonea into “early time signal™ and "late time signal", 
each of which characterize the signal by the time instant at 
(iimelmethe traditional model can be acceptable or not [Ref. 3]. 
This effort considers the development of a method that is 
capable of handling the “early time" whose time varying 
residues of exponentials do not become constant until the 


excitation of the incident field disappears. Such a signal 





can not be adequately or accurately modeled by the tradi- 
tional model. 

It is the intention of this work to investigate a method 
that can handle the "early time signal" and show "Robustness" 
under relatively heavy noise pollution as well as improving 
the accuracy of the poles we are interested in evaluating. 

A non-linear parameter optimization approach, through the 
modified least-squares minimization method, has been tried. 

in particular, we have been interested in evaluating its 
performance in improving the accuracy of the poles extracted. 
The non-linear parameter optimization approach using itera- 
tions has contribued to improvement of accuracy and showed 
the possibility of application to the "early time signal" as 
was the final goal. We define the "new signal model" as a 
Causal connection of the "early time signal" and "late time 
Signal" throughout this paper. The new signal model will be 
presented in Chapter II, with both the description of the 
transient response mechanism and the way in which the signal 
model is constructed through the transient response mechanism. 
Chapter III briefly describes the problems associated with the 
traditional model, and with many of the current extraction 
methods in existence. Chapter IV is devoted to the descrip- 
tion of the non-linear parameter optimization through a 
modified least-squares method and the way in which the algo- 
rithm is implemented. Chapter V presents test results of an 
attempt to use this approach for pole extraction and its 


performance characteristics. Chapter VI contains a summary 
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of the simulation results, the problems encountered, the 
potential problems expected, and the possibility of 
application of this method. 

Computations were performed using the IBM 3033 system at 
the Naval Postgraduate School. The testing of concept feasi- 
bility was the major impetus behind this effort, with concerns 
for practical implementation ina "real-time" environment left 
for subsequent study. Considerations to the processing speed 
will not be addressed in this thesis. All the calculations 


were done ina double precision environment. 


ll 





If. SCATTERING TRANSIENT RESPONSE SYSTEM 


A. TRANSIENT RESPONSE MECHANISM 

An electromagnetic wave incient upon a scattering body 
forces currents to be distributed such that Maxwell's equa- 
tions and the corresponding boundary conditions are satisfied. 
These induced currents produce a scattered field that can be 
defined at any coordinate in free space if the expressions for 
those currents are made correctly. Changes in the angle of 
incidence, the incident wave shape, the polarization mode, 
and the target geometry and composition affect the results of 
the evaluation of the scattered transient response. 

Basically, the singularities uniquely characterize a 
scatterer by their number and locations. The corresponding 
residues, which are interpreted as the weights of these poles, 
either can be viewed in the frequency domain or in the time 
domain. 

As shown in [Ref. 3], equation (2.1) represents the form 
of the scattered signal impulse response of the target, in 
digital form 


: Deel 
H (k) Hy (kT) + 2 Hy exp (kS_T) ( ) 


where Hg (t) =§ for t > Tf, T is the sampling interval and 
the infinite sum represents the natural resonance response 
with complex eae Oo . 3 jo in the left-half plane. The poles 


constitute the parameters for NCTR of an object. The process 


i 





of extracting those parameters from the impulse response 
first collects the data points within a finite time period 
(time window). 

The data collected is represented using equation (2.1) 
and equation (2.3) which partitions the signal into "early 


time" and “late time" components. 


00 M-1 
X(k) = } Hk) * J I(k) ee) 
k=0 k=0 
TZ-1 M-1 oo M-1 
X(k) = } H(k) * J) I(k) + } H(k) * JY Ik) (273) 
k=0 k=0 k=T9 k=9 
early time data late time data 


where I(k) is the incident wave at t = kT 


T is sampling interval in time. 


We can write the equation (2.3) in the form of equation (2.4) 
or (2.5). Both of these expressions display the essence of 


the new signal model which will be discussed in the next 


section. 
N 
mk) = (2A (k) exp (S_k) (2.4) 
where A, (x) = A, LOR Whee i 
N 
X(k) = E(k) + } A, exp(S_k) (225) 
gem IL 


where E(k) = 9 for k > T% 


and A, its constant for all k 
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T@ is the point at which the late time signal starts, so we 
consider that a time window might contain the data points 
ranging from the one before T# to that after TY. If we 
perform processing on a data window containing the data 
points which were from that early time region, we may lose 
some of the poles of natural resonance what are likely to be 
seen in the late time region. The data points after Kg 
presumably has all the information regarding poles, even 
though there are some extra poles that may be discriminated 
from the natural poles. These extra poles are introduced 
by the result of noise with relatively low SNR. Also, there 
is another risk to lose the natural poles having very high 
damping coefficients. There is no effective method yet found 
to compromise on having the point in which we can see all the 
poles with less extra poles under the higher SNR condition. 
Meanwhile, the discrimination of E(k) appears to us as 
the quantity to be handled to give improvement in accuracy 
of the natural poles. The next section describes the new 
signal model based on these aspects and the transient response 


mechanism. 


B. IMPULSE RESPONSE SIGNAL MODEL 
The late time impulse response can be expressed as a 
summation of exponentially damped sinusoids as in equation 


(ASIP A 


N 
mic) (2A exp(S_k), k=9,T,2T,..., (M-1)T (2.6) 
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where T = At, and M is the number of sampling points. 
The signal model expressed in equation (2.6) is valid after 
the incident wave has completely illuminated a target and 
no forced response remains. 

As the mechanism of a transient response system is strictly 
governed by the physics and orientation of the target and the 
attributes of the incident wave, we can write the equation 
(2.7) as an implicit form of equation (2.2), assuming we 
strike the target with a plane wave impulse. 


0 TZ-1 
X (k) = } H(k) - }  [ulk) - ul(k-T@) ] (27) 
b Kg 


@Gemaaning both equations (2.6) and (2.7), equation (2.9) would 


be sufficient if we assume N poles are present. 


wt.) = X_ (k) ++ Xk) 
00 TY-1 
X(k) = ) H(k) + }  [ulk) - u(k-T#) ] 
k=9 k=9 
~ A, exp(S_k), Ra Oe 2 eM) (229) 


We call the impulse response signal model, in the form of 
equation (2.9), an implicit form of the new signal model. 
A simpler form of equation (2.9) can be either of equations 


femreeyr or (2.11). 


N 
X(k) = E(k) + } A, exp(S_k) (2.19) 
n=1 
N 
X(k) = } A, (k)exp(S_k) (2s) 
n=l 


ike 


where A {k) = At ton ea oy 


In equation (2.1), the residues are time-varying until 
the time at which the augmented function E(k) in equation 
(2.19) becomes zero. Thereafter the A (x) = An (constant). 

To begin with equation (2.19), it is necessary to define 
a new independent variable as a parameter whose time behavior 
is not predictable along the positive time axis. It is con- 
venient to define this as a set of variables. By regarding 


them as independent, the equation (2.11) would be in the form 


of the equation (2.12) in the discrete digital data processing 


sense. 
N 
X(k) = E(k) + } A. exp(S_k) (Po NB} 
n=1 Nh nh 
TO-1 N 
Mie) € (k) + ) A exp (S_k) (Zo) 
n=l n=l 
assuming e,, (k) = Eo (k-n) 


The equation (2.14) is expressed in terms of a positive 
counting sequence in k. As in the equation, we might evaluate 
the equation (2.13) if we know the exact number of poles. But 
we know that it is impossible because we are in the region of 
the early time. So we write the equation (2.14) as a time 
shifted version such as the equation that has meaning from 
the observability viewpoint. 

TZ-1 


N 
X(k) = J} e.(k) + } A. exp(s_k) (2.14) 
n=l ™ n=1 ™ a 
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waeoe  K=-COtm, TYt+m-1l, ..., TH+1,TH,TG-1, ... , TE,G 
porem - 0 
Here, the examples of the new signal model are presented 
in Appendix C, in the form of a decomposed signal. This syn- 
thetically generated signal is written in terms of e,(k) 's 
and the sum of exponentials. In Chapter V, we present the 
decomposed signal from the synthetically generated data signal 


and that obtained by constructing the E(k) as the results 


of the non-linear parameter optimization processing. 


ee) NOISE 
White Gaussian noise is assumed throughout this work. By 
introducing the noise into the new impulse signal model, we 
call the following equation the modeling function of a tran- 
sient response system and the complete and general form of 
the new signal model. 
<0 


N 
X(k) = N(k) + } e (k) + } A, exp(S_k) (2 IL) 
n=l ” n=l 
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IIIT. PROBLEMS WITH THE TRADITIONAL METHODS 


In this chapter, a brief description of problems with 
the classic Prony's method, the Householder orthogonalization 
technique and the Eigenvalue method are provided. Each of 
these methods computes the location of poles in the s-plane 
but only the last two methods can estimate the number of poles 


prior to the discrimination of the poles. 


A. PRONY's ALGORITHM 
Equation (2.3) is interpreted as the following difference 
equation assigning the values of the real part of the signal at 


specified sampling points in its left side of the expression. 


; 
X(J) = A 
n=l m 
K(k) = AZ 
nal non (332) 
N 
X((M-1)k) = } AZ (M-1)k 


And a polynomial equation (3.2) that has the same roots 4, 


N 
) az = (2s 
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can be combined with equation (3.1) to yield the Prony's 
difference equation in the form of equation (3.3) 
N 


, Sees 0, tor k=0,1,2T, ... 
n=1 


Poe OE) ery (3.234) 


Then, we write equation (3.4) from equation (3.3) as 
N N-1 


) a X(n+k) = a_X(nt+tk) + a.x(n+k) = (cea 
n=l] ” wey = ‘ : | 


Using matrix notation, the equation (3.3) may be written as 


Selittron (3.5) or (3.6). 


N-lL 
L a X(n+k) = -a,)X(n+k) (3.5) 
Let ay = 1, then equation (3.5) becomes as 


ie N-1 ~ “Ny uo) 


where Xd is N by N circulant matrix of sampled data 
Aud is N by 1 Prony's coefficients matrix 


oy is N by 1 row matrix of sampled data set 


is + May) 


The Prony’s coefficient a's are to be calculated by 


making use of the characteristics of the circulant matrix X. 


= 
Semen Nef N-1' = N-1 “N cau 


As described in the above, this method can be applied after 
the number of poles are known to us, then the value of the 


s,'s are to be found from the equation (3.8). 
1g) 





eae LN(2_)/k, ois heme O20 (3.8) 


B. HOUSEHOLDER ORTHOGONALIZATION METHOD 
A general expression of equation (3.2) may be represented 


by equation (3.9) assuming the number of poles are not known. 
ag ~~ a,2 +... ta,Z = (3.9) 


where N' ais the unknown variable 
Then, the equation (3.9) is to be satisfied under the condi- 


tion that there are enough sampled data points m > 2N. 
agkXg + as} t+... Ft Qayikn = g (3.19) 


where the N by 1 sampled data matrix is defined as 


x, = PaGerr ee NCat2yy . 2 wy 6X {(M=N=141) J 


Application of the successive orthogonalization through the 


Gram-Schmit process would produce the orthogonal vector set. 


QO = Log, Se - e © 4 Ones 
n-1 
where o | = x, - iby <K_,0;> 0; Coe) 


If we set the Og toube lye enen Ene Weby 1 x slic igual GS) Ce(eeliglsy 


to be in the form of equation (3.11). 
PaO. ©) <X_, 0, > 9, | (Be 12) 


So that the equation (3.12) holds as X=0M. In the above, 


M is the 2-D matrix of multiplication factors whose diagonal 


20 





elements are all l's. Here, an orthogonal vector O71 makes 
the corresponding data set x, to be orthogonal to all the 
previous vectors Re - *i_2° If we have had all the x vectors 
during i-l times of successive orthogonalization process 
such that there is a non-zero orthogonal vector. We have to 
receive continuously the next data set Xeay and check whether 
the orthogonal vector in the next step is zero or under a given 
threshold. If it vanishes, we can say that we have (i-1l) 
poles. Our test against noise polluted data signal using the 
Householder method showed the inability of handling the early 


time signal as well as the data being heavily polluted. 


C. EIGENVALUE METHOD 
Again, from the equation (3.19), we can rewrite this 
equation as in the form of equation (3.14). 
ee = (354 ) 
q 


where se — [X91 * | . 8 « leenad 


The eigenvector can be derived from the equation (3.15) by 


doing some mathematical manipulations of equation (3.14). 


eee AS,” = Riga vane = g 


From this equation, we find the eigenvalues correspond to 
eigenvectors to see if there is any eigenvalue approaching 
zero. If there is one zero or under threshold, we also 


consider the number of poles as N'-l. Even both the 


Za 





Householder orthogonalization method and the Eigenvalue 
method can provide the means of calculating the number of 
poles, their algorithm (tending to fit the data points to 


the traditional signal model) showed the lack of generality. 
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IV. NON-LINEAR PARAMETER OPTIMIZATION APPROACH 


A. INTRODUCTION 

In 1963, Marquardt suggested an algorithm for the least- 
squares estimation of non-linear parameters [Ref. 6] when 
highly accurate parameter values are required. As we can 
see in the signal model in equation (2.14), that new signal 
model has multiple parameters that are functions of time. 
Here, we rewrite that equation again in the form of equation 


treks). 


TA-1 N 
X(k) = N(k) + } e.(k) + } A, exp(S_k) (4.1) 
n=k n=l 


where e, (kx) = §, for k >T@ 


We are not interested in the value of e, (kK) PESeIE, but 
include their estimations in order to contribute to the 
eeeuracy of s.'s and A'S: 

In other words, the extrapolation of the sum of exponen- 
tials, which is assumed to have all the poles, to points in 
the early time region may provide more accurate poles and 
residues simultaneously. 

The major advantage is the fact that we can make use of 
the high SNR data signal in the early time region with the 


basic information of poles that were derived in the late time 


region. 
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We will be processing on a data window whose vector 
length is to be increased one by one by moving the first 
element of that vector toward the time-origin point. The 
optimization process evaluates the normal equation to find 
the local optimized set and finds the global optimized set 
as its final goal. 

Let us define an ERROR function as in equation (4.3) 

N-1 


ERROR{A,,S_,E(k)} = 2 Penk) = X (ic) 


7 (Cee) 

where X'(k) is a measured data point 
To have minimized the least-squares error at every instant 
of measured time, the ERROR function has to be minimized by 
obtaining the global set of parameters being optimized. We 
also can write equation (4.3) in a more concrete manner as in 
equation (4.4). 

M-1 M-1 


y [x'(k) - x(k) ]*= J e(k)? (4.4) 
k=9 k=9 


ERROR {xX} 


where xX [E(k) ,A_,S_] 


Non-linear parameters such as e, (k) must be obtimized in 
a way such that the more accurate values of Sh have to be 
calculated as we increase the number of non-linear parameters 
e,; Users may define the numbers of parameter e,'s that 
are at least equal to or greater than the number of data 
points in the early time region. Now it is our task to find 
out the optimized value of Ansys and en through optimiza- 


tion processing, either in the global sense as it is represented 


in the equation (4.5) or using the normal equation (4.6). 
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GeenOoR(X) /oIX = G (4.5) 
de (kK) /o X(k) = g (4.6) 
where k=9,T,2T . , (M-1)K 
igi Beye , N 


With using the normal equation, as it is shown in the (4.6), 


a global set of non-linear set of parameters can be obtained 


through an iterative evaluation and fitting process. In our 


simulation work, the modified least-squares method was the 


basis of the equation (4.6). 


B. THE MODIFIED LEAST-SQUARES METHOD 


The error function is rewritten as in the expression of 


the equation (4.7). 
k 


ERROR{k,x} = X'(k) ~ X(k, x,) (4, 


where xX'(k) is the k-th sampled data 
kee x) 1s the model function 


xe is a vector containing the k-th sampled 


g 


Let us define xO to be an initial estimated value of x, 


a sequence of approximations to the optimized value is to 
generated by the equation (4.8). 
- m 
= - 4 
x = xX ee ea) J, ERROR(x ) ( 


where Jn is the numerical Jacobian matrix evaluated 


7) 


then 


be 


. 8) 


m is iteration number of successive optimization 


D. is a diagonal matrix equal to the diagonal of 


Zo 





an is a Marquardt parameter 


ime snumoer Of total iterations can be controlled by the 


threshold which is defined as in the equation (4.9). 


m+1 m (4.9) 
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V. TEST RESULTS AND PERFORMANCE EVALUATION 


A. INTRODUCTION 

In order to establish the ability of the program listed 
in Appendix B to extract the poles and residues correctly, 
three different simulated signals, each of which is polluted 
with infinite, 30 dB and 15 dB SNR, were created by the syn- 
thetic signal data generation routine. These 3 sets of signal 
data were chosen to span a range of possible situations. 

In the context of the transient signals and additive sta- 
tionary Gaussian noise being used, the definition of SNR is 
in terms of a ratio of energy quantities intergrated over the 


entire 20 nsec time window. 


TABLE I 


SIMULATED SIGNAL 1 


POLES 
RESIDUES (Nep. GHZ) 
AEE 1 fe game a ee ted 0 
sO coal. -1.9 - 11.9 


Simulated signal 2 consists of 2 sets of pairs of complex 
conjugate poles and residues, which is an extrapolation of the 
simulated signal 1, using its parameters. Simulated signal 3 
has 3 sets of pole-residue pairs, extrapolating from Signal 2. 

There are 5 options in the signal generation program in 
choosing a particular function of E(k). In this simulation 
work, a trapezoidal wave form of E(k) was used with TO=8 .42 


nsec. So that the unknown parameters that are to be augmented 
2 





at every processing step will be a maximum of 10 plus 4 times 
the number of poles at the final processing stage, when we 


have 512 data points within a 20 nsec time window. 


TABLE ITI 


SIMULATED SIGNAL 2 


POLES 
RESIDUES (Nep. GHZ) 
1.8 + il.g -1.9 + il.g 
1.9 - il.g -1.9 - i1.@ 
G.5 + ig.5 -2.8 + i2.9 
0.5 + if.5 -2.8 - 42.9 
TABLE III 
SIMULATED SIGNAL 3 : 
POLES 
RES IDUES (Nep. GHZ) 
1.9 + il.g -1.9 + il.g 
25 agen He Sa -1.8 - il .@ 
G.5 + if.5 -2.9 + i2.9 
0.5 - i9.5 ae ee PA | 
eee 5. + 16.25 -3.8 + i3.9 
9.25 - i9.25 -3.9 - 13.9 


The results of varying the number of additional data 
points from zero (all the data points are from late time 
region) to ten extra data points (10 extra points are from 
early time region are added to the points of late time 


region) are contained in Table IV through XII. These nine 
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@ases COrrespond to the 3 synthetic signals, each having 
3 additive noise levels. The 3 synthetic signals are plotted 
in Appendix C. In the next section, the accuracy of pole 
extractions, as indicated in the tables are displayed by way 
of graphical pole maps. In addition, reconstructed waveform 
obtained from the parameter extractions are compared graphi- 


cally to the original waveforms. 


Zo 
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PoratetiR OPTIMIZATION FOR SIGNAL 1 (NOISE Eee, i) 
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C. PERFORMANCE EVALUATION 

Unlike the traditional methods, this method revealed 
itself as one that can be applicable to a general scattering 
problem, including the early time data. A potential superiori- 
ty of this method over the traditional methods is found in its 
capability of handling the early time signal that has relative- 
ly large energy content and less sensitivity to the fixed noise 
level. 

Begarding the capability of handling the early time signal, 
the performance of this method seems to be very high compared 
with the traditional methods which do not model that signal 
generically. But, two points need to be considered when we 
ensure this high degree of performance. One important point 
is that we be able to estimate properly the number of poles to 
optimally ask for. This will depend upon the SNR of the data, 
increasing generally with increasing SNR. The actual physical 
scattering data has an infinite number of poles but, because 
the energy content of natural modes tends to drop off rapidly 
beyond the bandwith of the excitation signal the higher order 
poles are ignored, blanding with the noise present. Another 
point is that the initial data window has to contain only the 
late time signal data. 

As the results of the tables in the previous pages indi- 
cate, converging characteristics of the value of parameters 
is not blocked by a high level of noise pollution. It has 


been observed that the computed value of poles are likely to 


D2 





@emverge Closely to the true values that are listed in Tables 
feene=ough ITT either by choosing the number of poles as small 
or by increasing the size of the data window. Even though it 
imnane tO prove this optimality of this method (due to the 
non-linearity of the scattering transient response signal 
model), handling of the early time signal seems to be cast 
into a simpler problem space and that optimality may be ex- 
plained by observing the parameters at every processing step. 
Namely, the parameter estimation of this method is less 
influenced by the unreliable points than the more accurate 
ones. Figure 5.1 and 5.2 and the results in the previous 
section explain the reason why the above assertion is valid. 
We see that the estimated value of parameters through optimi- 
zation is not going to worse as we include the early time 
data points (these data points can be regarded as the unrelia- 
ble data when we extrapolate the sum of exponentials of the 
late time region into those of the early time region) into the 
data window what we see. But, as we are adding more data 
points of the late time region, improved values of parameters 
are obtained through processing, as in the results and Fig. 
5.2. Although this is true, we still see that the plot of 
E(k) in Fig. 5.2 does not provide sufficient data to confirm 
that an improvement of accuracy depends only on the number of 
late time data points, because those data points have a lower 
level of SNR locally. Now, let us get back to Tables IV, VII, 


and X. For the signal inputs whose SNRs are assumed to be 
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infinite (noise free), we have a set of parameters whose 
values are not converging to the true values as we add an 
increased number of unknowns. These results ensure us that 
the increased computer round-off errors obtained by more 
data processing (more parameters to extract) will produce 
divergence. Even initially with low SNR, the increased 
Signal energy in the additive early-time data may improve 
the results until, by adding enough extra data points and 
parameters, round-off error will eventually catch up with 
the most optimized one and beyond a specific point the 
results will begin to diverge. So that, in a general sense, 
we can say that the improvement in accuracy may be obtained 
either by the number of data points or by making use of the 
data points of high SNRs. If we consider a signal having 
high damping coefficients we can not use the signal data 
what are likely to be hidden by noise. The technique used 
here in processing the data proceeds by sliding the first 
data point along the time-axis towards the origin until the 
data window expands its range up to time-origin by adding 
the unknown parameters E(k) 's which is the same as the number 
of data points in the early time region. 

Figures 5.3 through 5.11 shows how rapidly the parameters 
we want to derive converge to the true values, within a con- 
vergence bounds. Here, all the plots were reconstructed from 


the results. 
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VI. SUMMARY 


A modified least-squares minimization method has been 
developed and tested against synthetically generated signal 
data using a non-linear parameter optimization algorithm. 

The effect of varying the number of data points either in 

eee region of early time or in the late time under various 
typical noise environments has been studied, defining the 
Semeeria for using this algorithm for the extraction of poles 
and residues under relatively heavy background noise, as well 
as for the time varying residue signals modeled by Morgan 
(Ref. 4]. 

In an attempt to use this method for signals with high 
damping coefficients, 3 synthetically generated signals were 
used as the inputs (Appendix C). Extraction of poles from 
the direct synthetic signal waveform was done through non- 
linear numerical evaluations, not by the Prony-based evalua- 
tions which were wsed in the traditional method. The 
parameter estimating algorithm for this non-linear signal 
model had worked successfully, providing users with optimized 
parameters that were converging to the true values, within 
specified convergence bounds. It was found that although we 
have the early time signal data having very high SNR, those 
data can not contribute to improving accuracy by simply 
increasing the number of extra data points and extra unknown 
parameters. 


36 





In order to get more accurate values of poles, we have 
to have an increased base of late time signal data, which 
have relatively low SNR. 

It has been shown that this method is "robust" even under 
heavy noise conditions. But three basic requirements have to 
be met for using this method as a general methodology; the 
optimal number of requested poles has to be known in advance 
to processing, reasonable initial estimation of parameters 
has to be made, and the transition time from early to late 


time signal models has to be known a priori. 
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APPENDIX A 
SYNTHETIC DATA GENERATION PROGRAM 
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